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SUMMARY 


The ACCESS-2 computer program was developed, based on 
successful experience with ACCESS- 1, to study the performance of 
approximation concepts when they are applied to problems of prac- 
tical size and complexity. Each modular process, such as equation 
solution, eigen- ana lysis or numerical optimization, is carried 
out m-core. As a consequence problem size is still limited by 
mam memory capacity, even though solvable problem size has been 
increased significantly, relative to ACCESS-1. Data transfer 
among various modules is carried out through sequential data files. 
Dynamic array allocation techniques are employed extensively to 
make the best use of available mam memory m each module. The 
introduction of dynamic array allocation leads to a more com- 
plicated program structure, however, the benefit of data com- 
paction is an essential ingredient of the quest for efficiency 
m the context of structural synthesis problems of practical size. 

A thermal load analysis capability is added with particular 
emphasis on the importance of thermal stress/strain considerations 
m design with fiber composite materials. This provides test 
cases for problems involving load vectors which depend on design 
variables. Frequency constraints are installed since they are 
important m their own right. In addition to being a good mea- 
sure of overall structural stiffness, natural frequencies repre- 
sent fundamental information that is essential to the formation 
of more advanced dynamic constraints (e.g. flutter) . Since direct 
application of ACCESS 1 type approximation (first order Taylor 
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series expansion with respect to the reciprocals of linked 
design variables) did not exhibit as good a performance as it 
did for static constraints, additional features for representing 
natural frequency constraints as first or second order Taylor 
series expansions with respect to regular linked variables were 
implemented. 

A constant strain triangular element with arbitrary ortho- 
tropic material properties was included m the element library 
to model laminated fiber composite material membrane structures. 
Membrane laminates are modeled by stacking elements of this type, 
with different material axis orientations, together. A thermal 
shear panel (TSP) element was introduced to take uniform soak 
temperature effects into account. A special shear panel element 
is needed because the symmetric shear panel (SSP) element for- 
mulation does not permit the neutral surface to change length, 
hence it cannot represent the thermal expansion associated with 
a uniform temperature change m the panel . The TSP element was 
especially devised ,for representing midplane symmetric wings 
subject to thermal soak temperature conditions m combination 
with aerodynamic and inertia loads. The temperature change 
effect is midplane symmetric and it is treated using the TSP 
element representation of the wing shear webs . On the other 
hand the mechanical loads on the wing induces midplane anti- 
symmetric response (bending and twisting) and it is treated us- 
ing regular SSP elements. The two results are then superimposed. 
Gradient information for the critical and near critical response 
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quantities is also computed by superposing independently 
computed symmetric and antisymmetric contributions. 

Optimization algorithms are completely independent of the 
other parts of the program and the current version includes the 
NEWSUMT optimizer. Only a few interface program statements 
would be required to add other optimizer options m parallel 
With NEWSUMT. 
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ACCESS - 2 


Approximation Concepts Code for 
"Efficient structural Synthesis 

User ' s Guide 


1. Introduction 

The development of the ACCESS-2 computer program was moti- 
vated by the successful demonstration, via the ACCESS-1 program, 
of the effectiveness of the approximation concepts m automated 
structural synthesis. The ACCESS-1 program was developed to 
test the performance of the coordinated use of approximation con- 
cepts on problems of relatively small scale, subject to simple 
static constraints. Furthermore, ACCESS-1 was designed as a 
base program that would lend itself to experimentation aimed at 
testing the effectiveness of new ideas and techniques for effi- 
cient structural synthesis. Therefore, the program and the data 
structures of ACCESS-1 were kept simple enough so that anyone 
with fundamental experience m computer programming would be able 
to understand the codes and modify them. As reported in Ref. 1 
efficiency, m terms of the number of finite element structural 
analyses needed to obtain near optimal designs, was improved 
significantly over previously reported capabilities having 
comparable generality. Furthermore, during the past two years 
at least four distinct research projects were carried out using 
ACCESS-1 as a base program and modifying it to test new ideas. 
Experience has shown that ACCESS-1 has been useful both m its 
own right and as a research tool designed to promote the use of 
approximation concepts in automated structural design. 
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However, many practical design problems are beyond the 
capacity of ACCESS-1 and consideration of more complicated con- 
straints than those treated m ACCESS-1 is often necessary and 
desirable. For example, m aeronautical applications, composite 
materials and thermal loads as well as stability and dynamic 
response constraints are extremely important. Indeed, there is 
frequently an even more compelling need for automated design 
capabilities when these rather complex materials and constraints 
are involved. The ACCESS-2 computer program has been developed 
m response to these needs and to build a body of experience that 
can be used to set effective guidelines for future development of 
large scale industrial application programs . 

The key ideas central to creation of the ACCESS-2 program 
are similar to those on which ACCESS-1 was built. Finite element 
analysis techniques and mathematical programming algorithms have 
been combined using a collection of innovative approximation con- 
cepts. Structures with prescribed configuration and given mat- 
erial properties are optimized so that their structural weight is 
minimized by modifying the sizing of finite elements; i.e. cross- 
sectional areas or thicknesses. (Optimum variation of configuratic 
and/or topological design features are still the subject of basic 
research investigations and no conclusive demonstrations of 
efficient design algorithms are available as yet) . 

The fundamental structure of the ACCESS -2 program is out- 
lined m Fig. 1. Upon activation, the preprocessor reads and 
prints out the input data m a readable format. The preprocessor 
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then computes all the ancxllary data that xs xndependent of 
changes xn the desxgn varxables and xt stores the results xn 
approprxate arrays as well as xn temporary external fxles (see 
Table 1) . When preprocessxng xs completed successfully, the 
desxgn process control (DPC) block xs actxvated and xt xnxtxalxzes 
the desxgn xteratxon process. At the outset the desxgn gxven xn 
the xnput data xs transferred to the approxxmate problem gener- 
ator (APG) , and thxs desxgn xs analyzed by the fxnxte element 
method. Constraxnt functxons are evaluated usxng the response 
quantxtxes obtaxned from the fxnite element analysxs and then 
the xnxtxal set of crxtxcal and potentxally crxtxcal constraxnts 
xs xdentxfxed and tagged. Explxcxt approxxmate expressxons for 
these tagged constraxnts are computed usxng the Taylor serxes 
expansxon wxth respect to approprxate xntermedxate desxgn varx- 
ables. Recxprocals of xndependent desxgn varxables are used as 
xntermedxate varxables throughout the program, except for an 
optxonal use of the xndependent desxgn varxables themselves when 
expandxng frequency constraxnts. In ACCESS-2, the objectxve 
functxon xs structural wexght and xt may be expressed exactly and 
explxcxtly xn terms of the xndependent desxgn varxables or thexr 
recxprocals . Thus , the APG block can generate an approxxmate 
problem statement of the form: 

Mxnxmxze W(X) 

£ - (x n ,x„ . . .x ) 

12 n 

Subject to 

H g (x) ^0 q = 1,2, ...Q 
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where W(X) and all H^{X) are explicit analytic functions of 
Note that the number of constraints Q for this approximate pro- 
blem is much smaller than that of the original structural design 
problem, because only the tagged constraints are included and 
all other constraints are temporarily ignored during a particular 
design stage. 

The data which define the approximate problem are sent back 
to DPC and subsequently given to the optimization algorithm block 
(OA) . The primary function of OA is to carry out numerical 
search process which will improve the design by operating on the 
current approximate problem statement. Since OA deals with 
problems that are stated m algebraically explicit form, it is 
not even aware that these problems are related to structural de- 
sign. Therefore, any established algorithm for inequality con- 
strained minimization of a function of many variables may be 
used. Ideally it would be desirable to select an OA from a 
collection of available options, taking into account the size and 
special functional characteristics of the current explicit 
approximate problem statement. However, m the current version 
of ACCESS-2 only the NEWSUMT optimization algorithm is available. 
This OA implements a sequence of unconstrained minimizations 
technique using a modified Newton's method and a quadratic extended 
penalty function feature to facilitate the unconstrained minimi- 
zations. One virtue of this interior penalty function type of 
OA (i.e. NEWSUMT) is that it can usually be controlled so as to 
provide an improved design that is also feasible with respect to 
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all of the constraint at each stage m the design process. 

After carrying out a numerical search with the approximate 
problem, the optimization algorithm (OA) block proposes an 

-y i 

improved design X to DPC. This step completes one stage 
of the design iteration procedure. 

In summary, one stage of iteration includes one finite 
element structural analysis, one constraint deletion process, 
one sensitivity evaluation for retained constraints, and one 
optimization of an approximate problem. Since the final 
design is subject to a detailed finite element analysis, the 
total number of finite element analyses equals the number of 
iteration stages plus one, which will be typically 10. The 
iterative design process is terminated when one of the speci- 
fied convergence criteria is satisfied. 

■+ » 

Strictly speaking each new design X proposed by the OA 
block to DPC is a better design than the original design X 
only with respect to the approximate problem statement. When 

• 4 -' 

the structural design corresponding to X is analyzed by means 
of the finite element method, it may turn out that some con- 
straints are violated. This situation may occur when the 
design changes m one stage exceed the applicable range of the 
approximate problem statement. It should be noted that the NEWSUMT 
optimizer is capable of locating feasible designs starting from 
an infeasible design. However, violation of constraints m the 
intermediate stages usually has a deleterious effect on the 
convergence characteristics and it is also inconvenient if 
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constraint violations are found at the final design. 
Constraint violation can usually be controlled or eliminated 
by appropriate use of the maximum step size parameter STEPMX 
and its dynamic modification feature via parameters STEPMX- 
multiplier and STEPMX-lower limit (see Appendix B) . 

All routines are written m standard FORTRAN XV language 
and they have been tested on: (a) the IBM 360/91 using the 
FORTRAN-H compiler at UCLA; and (b) the CDC 6600 at the NASA 
Langley Research Center. Implementation on other types of 
computers will be straightforward provided those computers 
have the required mam memory capacity. Except for the blank 
COMMON arrays, 380 -^qK and 280 ^qK bytes are required on IBM 
360/91 without and with program overlay, respectively. On a 
CDC 6600, the corresponding basic memory requirement is 
100gK words with overlay. 
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2 . Program Implementation 

ACCESS-2 computer program may be executed as a stand alone 
program. It consists of approximately 8500 FORTRAN statements. 

The program supplied upon request is a version operational either 
on IBM 360/370 or CDC 6600/7600 systems. Since it contains no 
machine dependent statements, it can be made operational on vari- 
ous computers, provided enough mam memory capacity and auxiliary 
data storage support are available. 

Auxiliary storage files are required as shown m Table 1. 
Files 10, 11, 12, 13, 14 and 15 are required for all problems. 

File 16 is required only when type 4 elements (TSP) are used m 
the structural model. Files 18, 19, 20, 21 and 22 are required 
only when second order expansions of frequency constraints are 
specified. 

The required size of blank common is very problem dependent: 
i.e. it depends on the structural analysis model (number of nodes, 
elements and load conditions) , the number of independent design 
variables, and the constraint types included. For certain pro- 
blems, it also depends on the initial design. Hence, it is rather 
difficult to give explicit formulas which estimate the size of 
blank common requirements . Table 2 gives actual blank common 
array size requirements for several example problems . 

Overlay or segmentation of the program can be designed 
easily by referring to Fig. 2. The simple 3 level overlay is 
adequate to solve most of the meaningful problems . If an oper- 
ating system allows more flexible overlay structure, it is pos- 
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sible to decrease the core requirement further. However, the 
net gam acquired by the elaborate overlay may not be signifi- 
cant, since most of the core is used for data and not for in- 
structions . 
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3. Structural Model and Input Data Preparation 

It is assumed that the reader is familiar with elastic 
structural analysis via the finite element displacement method, 
as well as with associated structural modelling techniques and 
typical data preparation procedures. Sufficient information for 
preparing the input data card images is given in Appendix B, 
therefore the explanations given in this section are limited to 
topics which require somewhat detailed technical discussion in 
order to avoid possible misunderstandings. 

3.1 System of Units 

Input data of the ACCESS-2 computer program may be prepared 

in any system of units as long as they are consistent. For 

example, if it is decided that the units for length and force are 

to be centimeters and Newtons, respectively, then the correspond- 

2 

ing units for pressure load or allowable stress must be N/cm . 
Note that the material constant specification calls for the spe- 
cific weight of the material, not its mass density. To be con- 
sistent lumped nodal mass should be given using weight rather 
than mass units. Example problems given in Appendix C are pre- 
sented both in the International System (IS) of Units and in the 
U.S. Customary (U.S.) units. Computer input data for examples 
are shown using numerical values associated with the U.S. units, 
simply because all the examples were originally presented in the 
literature using U.S. units. 

3.2 Node Numbering Scheme for the Effective Use of Memory 

The system stiffness and the mass matrices are stored in a 
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vector form within the skyline of the non-zero elements, i.e. 
there are no operations or no storage allocations with elements 
that remain zero during the solution. This scheme allows some- 
what more flexible node numbering arrangement than the ordinary 
band equation solver. It is better, however, to follow the same 
guidelines m preparing data as for a banded matrix solution 
scheme; i.e., differences among node numbers associated with an 
element must be kept as small as possible for all elements. 

3 . 3 Symmetric Wing Model 

If the webs of a midplane symmetric wing are modelled with 
SSP elements, only the upper (or lower) half of the wing is 
modelled. Assuming that the X-Y reference plane is the plane of 
symmetry, the X and Y displacement components and loading com- 
ponents are then anti -symmetric. Displacements and loadings m 
the Z direction are identical for both sides of the X-Y plane . 

For example, if a cantilever beam such as that shown m Fig. 3(a) 

is to be modelled using two SSP elements, then the simplified 
model should be that shown m Fig. 3 (b) . Note that only half of 

the load P need be applied to the node 3, since the other half 

is implicitly applied to the conjugate node 3’ (which does not 
exist explicitly m the model) . The SSP elements are always ver- 
tical to the X-Y plane of symmetry. 

The assumed' displacement function for SSP elements cannot 
accommodate uniform thermal expansion of each SSP element. If 
specified midplane symmetric temperature changes are specified 
for a midplane symmetric structural model, m which the vertical 
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webs are represented by SSP elements , ACCESS 2 branches and makes 
a separate calculation which adds m the midplane symmetric tem- 
perature change effects. This is accomplished by assembling 
equilibrium equations for the midplane symmetric structure with 
all of the SSP elements replaced by TSP elements while only con- 
sidering midplane symmetric temperature change loading. These 
equilibrium equations are solved for displacements u^ due only 
to midplane symmetric temperature changes . These thermally in- 
duced midplane symmetric displacements are superimposed on the 
previously computed midplane antisymmetric displacement state 
due to mechanical loads only. Treating the symmetric and anti- 
symmetric contributions separately reduces the number of dis- 
placement degrees of freedom that need to be considered m each 
of the two analyses and for thin wings it also tends to improve 
the accuracy of the analysis by avoiding the poor conditioning 
often associated with simultaneous treatment of bending and mem- 
brane response. The strain state is computed based on the total 
displacement/ and the stress state is computed by transforming the 
strain state using the stress-strain relationships. 

3.4 Design Variable Linking 

The general concept of design variable linking is discussed 
m Sec. 2.3.1 of Ref. 1. In the ACCESS-2 computer program, if 
the sizes of some group of finite elements of the same type are 
controlled by a single design variable, these elements are said 
to belong to the same design variable linking group. The sizes 
of elements m a design variable linking group are modified m 
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proportion to the initial sizes given in the input data. 

Design variable linking groups are also used to define "re- 
gions" for the regionalization of stress constraints. The general 
idea of regionalization is described m Sec. 2.4.1 of Ref. 1. 
Elements which belong to the same design variable linking group 
form a region and only one stress constraint per load condition 
(the most critical) is considered for each group m any stage of 
the iterative design procedure. Selection of the critical stress 
constraints within a region is not rigidly fixed, but dynamically 
updated at the beginning of each stage. If the location of the 
critical stress constraints shifts frequently within a region 
from stage to stage the iteration process may be unstable, al- 
though this type of instability was not observed in solving any 
of the problems given m Ref. 1. However, if the user desires to 
remove the regionalization of stress constraints, it is only 
necessary to specify IGLINK = -200. 

3 . 5 Failure Criteria for CSTOR elements 

The CSTOR element is implemented to model structures made 
with orthotropic materials including multi-layered fiber composite 
laminates. While strength failure criteria for isotropic metal 
alloy materials are imposed using the von Mises combined effective 
stress, strength failure criteria for CSTOR elements are selected 
from 3 available options . 

They are: 
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A. Maximum strain criteria 


S L * e L - “l 4T £ 

^ e T - a T AT < ej 


T LT ' “ Y LT 


B. Stress interaction formulas 



£ 1 


C. Tsai-Azzi Criterion 




< 1 


where 


'LT 




T 


LT 


-c 

£ L 


longitudinal strain 
transverse strain 
shear strain 
longitudinal stress 
transverse stress 
shear stress 

allowable longitudinal compressive strain 
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: allowable longitudinal tensile strain 

Y lt : allowable shear strain 

F ,F : allowable longitudinal and transverse stresses 

Li i 

F Tfn : allowable shear stress 

JLjL 

v TL : Poisson's ratio relating to contraction m 

1 " the longitudinal direction due to extension 

m the m-plane transverse direction 

v LT : Poisson's ratio relating to contraction 

in the m-plane transverse direction due 
to extension m the longitudinal direction 

Among the three alternative strength criteria, the maximum 
strain criterion is the most conservative while the stress inter- 
action formulas are usually the least conservative. 

3 . 6 Computation of Constraints 

All constraints*, except the side constraints, are normalized 
so that potentially critical constraint functions m the feasible 
region assume values between 0.0 and 1.0. Constraint functions 
are defined as follows: 
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Slope (Relative Displacement) Constraints 
Slope 

S (U) - (S 0 -S.)/d > 0 

s lu) 

Relative Displacement 

r (U) - (6 2 -6 1 ) > 0 
r 

where d is the projection of the distance between the two points 

ir 

on a plane normal to the displacement components <5^ and 5 ^. 


Stress (Strain) Constraints 


a (U) - a 

_ CO) 


i 0 


a - a (L > 

TuT - 


^ o 


For strain constraints, see 3.5, 


Frequency Constraints 


(U) 2 

to' ' ~ to 


to 


(U) 


> 0 


u 2 - » (L) 


to 


(L) 
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4. 


Restrictions and Limitations 


The amount of main memory storage required for solution of 
a particular problem depends upon many factors, including the 
number of nodes, the number of elements, the number of design 
variables, the element types used, the hinds of constraints im- 
posed, and even the initial design employed etc. For static 
problems, it is necessary to retain two system stiffness matrices 

and the load vectors in core. For dynamic problems, three system 

* 

matrices must be retained m core. If a problem involves dynamic 
constraints and thermal shear panel elements, four system matrices 
must be in core simultaneously. Also a complete approximate pro- 
blem statement (all retained constraint values and all the corres- 
ponding derivative components) must be m core for the OA block. 

It is difficult to estimate the array size required for a system 
stiffness matrix m advance. Only the nonzero skyline of an upper 
half matrix is stored, hence the memory requirement depends on 
the node numbering scheme. For medium size problems (300 -600 
DOF) , the density of nonzero elements in the matrix is usually 
20-50% and a first approximation can be made by estimating the 
density based on observation of the finite element model. The mam 
memory storage required for the integer portion of the blank 
COMMON is usually less than 10,000 words, but the real variable 
portion is very dependent on the nature of the problem. For pro- 
blems m which the number of constraints retained tends to be 
larger and in which there are many independent design variables 
(e.g. structures involving laminated fiber composite skins) the 
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constraint derivative array size [i.e. (Number of Design Vari- 
ables) x (Number of Constraints Retained) ] may limit the problem 
size, since this large array must be m core m addition to the 
instructions and local variables. 

When first order approximations are used frequency con- 
straints can be imposed on any subset of frequencies within the 
lowest NFREQ frequencies. If second order approximations are 
employed all frequencies m the lowest NFREQ frequencies must be 
bounded. 

Capabilities for aeroelastic constraints are not available 
m this version, therefore NMODE must be zero and the flight 
condition specification flags must all be zero. 

All input data are read m with fixed format, hence column 
positions of the punched data are of critical importance. Es- 
pecially note that all blank columns are regarded as zeroes for 
numerical inputs. 
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Table 1. Temporary Files 


File Name 


Contents 


10 

11 

12 

13 

14 

15 

16 

17 

18 

19 

20 

21 

22 


Stiffness matrix components associated with 
unit values of independent design variables 
and load vector components which are mdependen 
of design variables 

Mass matrix components associated with unit 
values of independent design variables 

Load vector components due to thermal loads 
and dependent on independent design variables 

Load vector components due to inertia loads 
and dependent on independent design variables 

Constraint gradients 

Input data and a part of preprocessor output 

Thermal shear panel stiffness matrix components 
Used only when IETP(6)^0 

Not used 

Eigenvector sensitivity vectors, if computed 

Mass matrix post-multiplied by eigenvectors. 
Required only when second order expansion of 
frequency constraints is used. 

Original system stiffness and mass matrices. 
Required only when second order expansion of 
frequency constraints is used. 

Modified [K-X-lM] m the eigenvector sensitivity 
computation. Stored m decomposed form. 

- L L ’ ((1=1 ' NE1G) ' b=1 ’ B) 

Required only when second order expansion of 
frequency constraints is used 
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Problems 


Wing Carry-Through 
Truss Model 
(static) 

Delta W] ng (Metal) 
(Static & Dynamic) 

Delta Wing-Composite 
(Static & Dynamic) 

Delta Wing-Composite 
(Static, thermal 
and Dynamic) 

Elevon-Composite 
(Static, thermal 
and Dynamic) 


Elements 

Total No. 

Free Displ. 

No. of 

Total No. of 

Required 


of Elements 

d.o.f -s 

Design Variables 

Constraints 

Real Array 

Integer Array 

TRUSS 

63 

42 

* 

63 

319 

11415 

1896 

CST 

63 

105 

28 

164 

7400 

2033 

SSP 

70 






CSTOR 

SSP 

252 

70 

105 

60 

2725 

12750 

6747 

CSTOR 

252 






SSP 

70 

105 

60 

2725 

16960 

8264 

TSP 

70 






TRUSS 

45 






CSTOR 

SSP 

192 

45 

109 

68 

3354 

47334 

10680 

TSP 

45 











APPENDIX A 


ELEMENT LIBRARY 


Currently, 6 element types are available: they are TRUSS, 
CSTIS , CSTOR, SSP , PSP and TSP. Basic characteristics of these 
elements are given in the sequel. 

1. Type 1 - TRUSS ; Pin jointed bar element of uniform cross 

section 



Strain-Displacement Relation (local coordinate) 


e 



1 ] 



(A-l) 
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Stress Strain Relation (local coordinate) 


a = E e - EaAT (A-2) 

where a : thermal expansion coefficient 

AT : average temperature change 


Force Displacement Relation (local coordinate) 


E A 
L 


1 -1 
-1 1 


u_ 


u 


QJ 


EaATA 




(A- 3) 


where F , F are externally applied force at P and Q 
nodes, respectively. 


Force Displacement Relation (reference coordinates) 
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Consistent Mass Matrix (reference coordinates) 


[M] = 


pAL 

6 


0 0 
2 0 
2 

Sym 


1 0 
0 1 
0 0 
2 0 
2 


0 

0 

1 

0 

0 

2 


(A-5) 


where p : density 
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Stress-Strain Relation (local coordinate) 



Stress-Displacement Relation (local coordinate) 



Local-Reference Displacement Relation 
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: unit vector parallel to the x-axis 


where 


1 


x 


1 


: unit vector parallel to the y-axis 


Up,UQ/U^ : displacement vectors of P,Q,R nodes. 


[T] = 


£„ m n 
xxx 


£ m n 

y y Y 


£ m n 
xxx 


£ m n 

y y y 


£ m n 

xxx 

£ m n 

y y y 


(A-10) 


Stiffness Matrix (local coordinate system) 


K = K + K 
n s 


(A-ll) 


where 


K 


n 


Et 

4A(1-v 2 ) 


(s-b)^ -v(s-b)h 

- (s-b) s 

v (s-b) h 

(s-b)b 

0 

h 2 

vhs 

-h 2 

-vhb 

0 


2 

s 

-vhs 

-bs 

0 

Symm. 


h 2 

vbh 

0 
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where p : density 
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hs 


0 -bh 


(s-b)s 0 (s-b)b 
-hs 0 bh 

s 2 0 -bs 

0 0 
b 2 

coordinates) 

= f (A-12) 



(A-13) 
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Stress-Strain Relation (material axis) 



v 

1— v 


tl e l 

LT V TL 


E. 


1 -' , LT V TL 



W'WV 

1 - V LT V TL 




E T( a T +V LT a L) 

1_V LT V TL 






(A-15) 


Strain Transformation Law (material-local) 



l 

l 
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2 

Ly 

2 

Ty 

Ly^Tx 

[T] 


^Lx^Ly 

ft Tx*Ty 

^Ly^Tx 


+ £ 


Ty^Lx 



where l. - 
Lx 

-T 
e L * 

e = cos 0 , 

X 

^Tx 

_ AD 
e T 

- = -sin 

Ly 

AT 
S L ‘ 

e" = sm0 , 
Y 

& Ty 

AT 

= e rr> 

T 

. e“ = cos 6 
Y 


Note: the direction of is chosen so that 

(e^ x ^ 2 ) • C? x x IT ) > 0. 


(A— 16 ) 


Stress-Displacement Relation 



[D] [T] [B] u - AT • li 


(A-17) 
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Local-Reference Displacement Relation 
same as type 2 

Stiffness Matrix (local coordinate) 



v lt e t v tl e l 


p = smB 
U = cos 8 
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(s-b) D x 


+ 2h(s-b)D 0 +ht D 


h (s-b) D-^ 

h 2 D 1 -2h(s-b)D 2 

+ [h 2 - (s-b) 2 ]D 2 -h(s-b)D 3 

+(s-b) 2 D 3 

-s (s-b) 

-hsD 1 [ (s-b) s-h 2 ]D 2 

-(2s-b)hD 2 -h 2 D 3 

+ (s-b)h D 3 

-h(s -b)D^ 

-h 2 D 1 +h (2s-b)D 2 

+[(s-b)s-h 2 3D 2 +sh D 3 

-s (s-b) D 3 

b(s-b)D 1 +bh d 2 

bhD^-b (s-b) D 2 


bCs-bJD^-bh 


-bhD 2 +b(s-b)D 3 


where D, = 


4p 2 U 2 

2pU (p 2 -y 2 ) 

2 2 2 
(p -y 2 ) 










Equilibrium Equation (local coordinate) 


K u + 


s = t 


h = t AT V 


- (b-s) (p 2 h 1 +y 2 h 2 ) + 2hpy(h 1 -h 2 ) 
-h(y 2 h 1 + p 2 h 2 ) + 2 (b-s) py (h^-h 2 ) 
-s(p 2 h 1 + y 2 h 2 ) -2hpy (h^-h 2 ) 
h(y 2 h 1 + p 2 h 2 ) + 2spy(h 1 "h 2 ) 
b(p 2 h 1 + y 2 h 2 ) 

-2bpy (h-^-h 2 ) 


> 


(A-21) 


where 


h - E L (a L + V TLV 

1 ^” V LT V TL 


h 2 


V a T + V LT a L ) 
1 “ V LT V TL 


Consistent Mass Matrix (reference coordinate) 
same as type 2 
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4. Type 4 - SSP: Symmetric shear panel element with uniform 

thickness and isotropic material 

This is a special element used to model relatively thin 
symmetric structures such as idealized supersonic lifting 
surfaces. Theoretical discussion is given m Ref. 1. It 
is assumed that this element models the upper (or lower but 
not both) half of the symmetric structure and the element 
plane of symmetry coincides with the X-Y plane. It is 
further assumed that all SSP elements are placed vertically 
with respect to the X— Y reference coordinate plane. 



1. There are only two nodes per element. 

2. The line of intersection with the XY plane does not 
move in the XY plane. It can only move vertically. 
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3 


. If the heights PP' and QQ' are different, the average 
(PP' + QQ')/2 is considered as the height of the 
element, i.e. b. 

4 . No thermal load can be considered m this element 


Strain Displacement Relation {local coordinate) 



e = [B]u 
wherein n = Y/b 


(A-22 ) 


Stress-Strain Relation (local coordinate) 


x 


E 


1-v 


xy / 


1 

V 

0 

e \ 

V 

1 

0 

/ X | 

0 

0 

1-v 

< Y / 



2 

> 

Vxy/ 


(A-23) 



Stress Displacement Relation (local coordinate) 
s 



_ 2r\ 
a 

0 

2n 

a 

0 

0 

0 

0 

0 

1 

1 

1 

1 

2 (1+v) b 

2 (1+v) a 

2 (1+v) b 

2 (1+v) a 


\ 


/ 



(A-24) 
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SXZTY 


Local to Reference Displacement Transformation 

\ 


u T 


Pi- 


ll 


£ 


x 

0 


m 0 
x 

0 ± 1 * 


0 


£ m 
x x 


Q 


Q' 


* (-) sign if Zp < 0 and < 0 


0 

± 1 * 


U T 


V„ 


W T 


U 


V 


\ W Q 


(A-25) 


where £ and m are components of a unit vector e along 
x x ^ x 3 

the local x axis. 


Stiffness Matrix (local coordinate) 


K = 


Et 

12 (1+v) 


F+3a -3 -F+3a 3 


3/a -3 -3/a 


Symm. 


F+3a 


3/a 


(A— 26 ) 


where a = 


a 

b 


F = 


2 (1+v) 


a 
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Consistent Mass Matrix (local coordinate) 




| + G 





0 


M = P^t 
2 



1 

6 


0 


1 

9 


0 


0 


1 

18 


(A-27) 


1 

9 


+ G 



1 

3 


Symm. 


1 

3 


0 

1 

9 


where p = 


density 
2 


r _ _v_ v 

^ “ 30 18 


30a 


H 


a 


+ 


1 

a 


It may look strange that the mass matrix depends upon 
Poisson's ratio v through G. This is due to the fact that 
the assumed displacement field is derived based on assumed 
stress field. (see Refs. 1 and 3) 
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5 


Type 5 - PSP : Pure symmetric shear panel element with uniform 
thickness and isotropic material 


This element is identical to a type 4 (SSP) element, 

except for a minor change m the assumed displacement state 

so that the stress state of the element is pure shear: i.e. 

o = a = 0 . This implies that e = e. = 0 . 
x y ^ x y 


Strain Displacement Relation (local coordinate) 
Y xy =[e,- li.lj {v V V *(jf 

Stress-Strain Relation (local coordinate) 


E 


T xy 2 (1+v) ^ 


xy 


Stress-Displacement Relation 
E 


xy 


2 (1+v) 


b' a' b' a] { U P f V P' U Q' v q|. 


Local to Reference Displacement Transformation 


(A-28) 


(A-29) 


(A-30) 


same as type 4. 

Stiffness Matrix (local coordinate) 
same as type 4 except F s 0 . 


Mass Matrix 

Assumed to be the same as type 4. 
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6. Type 6 - TSP : Thermal symmetric shear panel element with 

uniform thickness and isotropic material 

Since SSP and PSP cannot be used for problems involving 
thermal loads, this special element is added to the ACCESS-2 
element library. The TSP element is designed to be used under 
steady thermal soak load conditions such that the temperature 
change in each TSP element is uniform and therefore symmetric 
with respect to the X-Y plane. 

If the structure is subject to both mechanical and thermal 
loads, two structural models must be created and analyzed separ- 
ately. One model is to use SSP elements to model shear panels 
and it is subject to only mechanical loads. The other model uses 
TSP elements to model the shear panels and it is subject to only 
thermal soak loads. These two models are created automatically, 

if the user specifies both SSP and TSP elements. Displacement 

♦ 

and stress states of the structure subject to both thermal and 
mechanical loads are generated by superimposing the results ob- 
tained from the two separate models. 

Theoretically, it is also possible to consider the PSP - TSP 
element combination, but this is not implemented m the current 
version of ACCESS-2. 

Note that the TSP option requires a significant amount of 
core memory and CPU time, since two system stiffness matrices 
are stored and decomposed. Sensitivity analyses of the responses 
must be carried out separately and superimposed afterwards. 
Therefore, analyses effort is nearly doubled when thermal effects 
need to be considered. 
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ORIGINAL PAGE IS 

OF POOR QUALITY 


Strain-Displacement Relation (local coordinates) 




Stress-Displaceitient Relation 



Local to reference displacement transformation 


same as type 4 
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Stiffness Matrix (local coordinate 



where a = §• 
b 

Force Displacement Relation (local coordinate) 



(A-34 ) 


(A-35) 


Consistent Mass Matrix - 

/ 

Assumed to be the same as type 4 

Note: As shown m the stress-displacement relation, stress 

distribution is linear with respect to x. In order 
to simplify the problem, an approximate stress dis- 
placement relation is used in computing stress and 
stress sensitivity. 
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ORIGINAL, PAGE IS 
OF POOR ' QUALITY 
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APPENDIX B 


INPUT DATA DESCRIPTION 

I. Job description and heading cards (II, 79A1) 

The first column is used as follows 

0 or blank : ordinary heading cards , whose contents m 

columns 2-80 will be printed on the first 
page of the output. 

1 : indicates that this is the last heading card 

and input data cards follow. 

2 : request for immediate normal termination of 

this job. 

Any number of cards may be used to describe or to comment 
the job. No.te that the last heading card must have "1" 
punched m the first column. Without this, all of the 
data may be regarded as heading cards . 

II. Primary control cards 

Card 1 (715) 

IOPT : 1 = Input data check only 

2 = Structural analysis only 

3 = Structural analysis and constraint 

function evaluation 

4 = optimization by the NEWSUMT optimizer 

5 

* ’ Not used yet 

IPRINT ; Printout control parameter except for output 

from each optimizer. 

2 = Standard output 
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IGLINK : 0 = standard execution 

-200 = removal of stress constraint regionali- 
zation 

-300 = removal of stress constraint regionali- 
zation for fixed size elements only 

IANALY { 1 ) -IANALY ( 4 ) 

: Not required to specify (leave as blanks) if 
IOPT z 4. Must specify if IOPT r 3. 

IANALY (1) : 1 = Compute displacement 

0 = Skip displacement calculation 
IANALY (2) : 1 - Compute stress/strain for all elements 

0 = Skip stress/strain calculation 
IANALY (3) : 1 = Compute eigenanalysis 

0 = Skip eigenanalysis 
IANALY (4) : 0 always 

Card 2 (1015) 

IN : Total number of nodes 

IBN : Number of boundary nodes 

INL : Number of load conditions 

IMATIS : Number of isotropic materials 

IMATOR : Number of orthotropic materials 

INITVG : Number of initial value groups for design varia 

ILOWBG : Number of minimum size groups for design variat 

IUPPBG : Number of maximum size groups for design vanat 

ITHLDG : Number of thermal load groups 

IPRLDG • : Number of pressure load groups 
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Card 3 


(1015) 


.IETP(i) , i = 1,2, - . .10 

Number of elements in the i-th element type 
i = 1 TRUSS 

2 CSTIS CST isotropic 

3 CSTOR CST orthotropic 

4 SSP symmetric shear panel 

5 PSP pure shear panel 

6 TSP thermal shear panel 

III. Node Coordinates (15, 5X, 3E10.4) 

IN cards are required to specify the node coordinates of node 
numbers 1 through IN. The order of the cards may be random. 


n i 

: Number of the i- 

-th node 


xn l 

: X 

coordinate 

of 

the 

node 

n 

Yn 

: Y 

coordinate 

of 

the 

node 

n 

Zn. 

: Z 

coordinate 

of 

the 

node 

n 


IV. Boundary Conditions (415, 3E10.4) 

If all 3 degrees of freedom' associated with a node are free, 
the node is not a boundary node . Otherwise it is a boundary 
node and for each boundary node, a card is required. 


fc> n i 


i-th boundary node 

number 

IBX bn 

n 


constraint code: 0 

— free 

J. 

IBY, 

bn 


- 1 

= fixed 

IBZ bn 

l 

/ 

2 

* = prescribed nonzero 
displacement 


* not available in the current version 
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prescribed nonzero displacements 
required only for constrain code = 2 


Element Data 

If IETP(i) = 0 for the l-th element type, no data is required. 
For each element type with IEPT(i)^ 0, IETP(i)+l cards are 
required. 


Card 1 


element type number (15) 


Card 2-IETP(i)+l s element information (1115) 

M : element number 

NP : node number corresponding to the internal 

node number P 


node number Q 
node number R 


node number S 


linking group number, = 0 for the fixed size 


LBGN 


UBGN 


MTLGN 


elements 

initial value group number 
lower bound group number 
upper bound group member 
material group number 
> 0 for isotropic materials : 1,2,... 

< 0 for orthotropic materials: -1,-2,... 
side constraint code 

-1 : element size restricted by the lower 


bound only 
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0 


: non negativity constrain only 

1 : element size restricted by the 

upper bound only 

2 : element size restricted both by lower 

and upper bounds 

Comments 

1. Elements must be numbered starting from 1 through 
IETP(i) for each element type. For example, if a 
structure is modeled using 100 TRUSS elements and 

300 CST elements, TRUSS element numbers are 1,2, 3,... 100 
and CST element numbers are 1,2, 3... 300. Within an 
element type, order of element data cards may be random. 

2. NR and/or NS are not required for element types with 
only 2 or 3 nodes per element. 

3. LGN, linked group number starts from 1 for each element 
type. For example, if a structure is modeled with 

100 TRUSS and 300 CST elements, with 10 and 30 design 
variables allocated to TRUSS and CST, respectively, then 
the linked group number for TRUSS runs from 1 through 10 
and that for CST ranges from 1 through 30. 

VI. initial Values (7E10.4) 

INITVG real numbers must be given. If INITVG > 7, two or 
more cards are required. The first value of the first card 
indicates the initial value for the group number 1, and so on. 

VII. Lower Bound Values (7E10.4) 

Minimum gauge values. ILOWBG real numbers must be given. 
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If ILOWBG > 7, two or more cards are required. If 
ILOWBG = 0, no card is required. 

VIII. Upper Bound Values (7E10.4) 

Maximum gauge values. IUPPBG real numbers must be given. 

If IUPPBG > 7, two or more cards are required. If 
IUPPBG = 0, no card is required. 

IX. Isotropic Material Data (6E10.4) 

IMATIS cards are required and on each card the following 
6 real numbers must be given. 

E : Elastic modulus 

v : Poisson's ratio 

Y : Specific weight 

a : Thermal expansion coefficient 

cr TD : Allowable compression stress 

.LJd 

a rT _ : Allowable tensile stress 

UB 

X. Orthotropic Material Data (7E10.4/7E10.4/6E10.4) 

IMATOR X 3 cards are required, i.e. for each material group 
3 cards are required, containing the following data. 

Card 1 

E : Longitudinal elastic modulus 

Xj 

E t : Transverse elastic modulus 

G t _ : Shear modulus 

-LiX 

v Tm : Longitudinal Poisson's ratio 

XjX 

y : Specific weight 

a T : Longitudinal thermal expansion coefficient 

Xj 

a T : Transverse thermal expansion coefficient 
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XI. 


Card 2 




in 


L 

t 

L 

C 

L 

t 

T 

c 

T 

LT 

t 

L 

c 

L 

t 

T 

c 

T 


LT 
Comments 


Direction cosines of the longitudinal 
axis with respect to system reference 
coordinates 

Tensile allowable longitudinal strain 
Compressive allowable longitudinal strain 
Tensile allowable transverse strain 
Compressive allowable transverse strain 
Allowable shear strain 
Tensile allowable longitudinal stress 
Compressive allowable longitudinal stress 
Tensile allowable transverse stress 
Compressive allowable transverse stress 
Shear allowable stress 


1. The transverse Poisson's ratio v__ is internally 

XXj 

computed using the relation v rPT E T = v rm E„ 

2. Depending upon the failure criteria applied to the 
specific material, either strain allowables or stress 
allowables are left unspecified. Failure criteria 
options will be specified later m the category XX. 

Lumped Nodal Loads 

Two card groups are required to specify lumped nodal loads 
applied to the structure. 
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Card Group 1 (1415) 

Number of nodes subject to lumped nodal loads for each 
load conditions. 1NL integer numbers must be given. 

Card Group 2 (15, 5X, 3E10.4) 

For each load condition, the specified number (by the 
group 1 cards) of cards must be given to identify the 
node numbers and associated load components m the 
reference coordinate system. 

XII. Pressure Load Data 

No card is required if IPRLDG =0. If IPRLDG > 0, the 
following 5 groups of cards must be given. 

Card Group 1 (1015) 

Number of elements subject to pressure load for each 
element type. (Presently, only CSTIS and CSTOR elements 
can be subject to pressure loads) . 

Card Group 2 (1415) 

Pressure load ON-OFF flag for each load condition, 
k 

ONOFF = 0 No pressure load for load condition k 
- 1 Pressure load should be considered for 
the k-th load condition. 

Card Group 3 (1415) 

Element numbers subject to pressure loads for all member 

MTYP 

types corresponding to NEPL 7 * 0. For each element 

type, the first element number subject to pressure load 
must be punched in columns 1-5; namely the group 3 cards 
should be subgrouped for different element types. 
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Card Group 4 


(1415) 


For each load condition corresponding to a load condition 
with ONOFF = 1/ an identical amount of data similar to 
that specified m the card group 3 must be given. Those 
' numbers designate the pressure magnitude group numbers , 
which are the pointers to the pressure magnitude applied 
to the corresponding element type and element number. 

This set of cards should be given for all load conditions 
with ONOFF^ = 1. 

Card Group 5 (7E10.4) 

Pressure load magnitude for each pressure load group 
must be given. IPRLDG real numbers are required. 

Comments : 

1. The direction of the pressure force is determined 

by the node numbering scheme of the triangular element 
and also by the sign of the pressure load magnitude 
specified m the card group 5. When the P, Q and R 
nodes of the triangle are in counter clockwise order 
and the corresponding pressure magnitude has a posi- 
tive sign, positive pressure is applied to the surface 
of the triangular region. 




Fig. B-l Pressure Load Sign Convention 
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2. Pressure applied on a single triangular surface must 
be uniform; no variation of pressure over an elements 
surface can be represented. 

XIII. Inertia Load Data 

Self-weight m a gravitational field or uniform translational 
acceleration will be accounted for by specifying this set of 
data. Note that rotational inertia loads cannot be considerei 
Two groups of cards are required. 

Card Group 1 (1415) 

INERTL k : k * 1 , 2 , . . . INL 

Intena load ON-OFF for each load condition. 

1 : Inertia load exists 

2 : No inertia load for the load condition 

Card Group 2 (4E10.4) 

For each load condition with INERTL ^ 0 , one card will 
be required. 

Tr 

ACC : Magnitude of acceleration in units of the 

standard earth gravitational field (i.e. 4g) 

X 

Direction cosine components of the acceleratio 
vector in the reference coordinate system. 

Z sj 

XIV. Thermal Load Data 

no card is required if ithldg =0. if ITHLDG > 0, the 
following 5 groups of cards must be given. 

Card Group 1 (1015) 

Number of elements subject to thermal load for each 
element type. 
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Card Group 2 (1415) 


Thermal load ON-OPF flag for each load condition. 

V 

ON-OFF =0 No thermal load for load condition k 

1 Thermal load should be considered for the 
k-th load condition. 

Card Group 3 (1415) 

Element numbers subject to thermal loads for all member 

MTYP 

types corresponding to NETH ^ 0 . For each element 
type, the first element number subject to thermal load 
must be punched m columns 1-5; namely the group 3 cards 
should be subgrouped for different element types. 

Card Group 4 (1415) 

For each load condition corresponding to the load condition 
with ON-OFF = 1, an identical amount of data similar to 
that specified m the card group 3 must be given. Those 
numbers designate the temperature magnitude applied to 
the corresponding element type and element number. This 
set of cards should be given for all load conditions with 
ONOFF k = 1. 

Card Group 5 (7E10.4) 

Temperature change for each thermal load group must be 
given. ITHLDG real numbers are required. 

Comments ; 

1. Each element is considered to have uniform temperature. 

2. Temperature change should be computed with respect to 
an appropriate uniform reference temperature. Note 
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that if all elements are made of the same material 
and assume the same temperature, then thermal stress 
are not induced. 

XV. Flight Condition Data 

This block of data will be reserved for future development 
of ACCESS-2 program which may include flutter constrains. 
Specify 0 for all IFLIGH k , k = l,2,.v.INL. 

No additional cards are required. 

XVI . Lumped Nodal Mass Data 

Card 1 (15) 

NMASS : Number of lumped nodal masses 

Card 2~ (NMASS+1) (15, 5X, E10.4) 

: Node number to which the mass is attached. 

: Weight of the mass 

Comment : 

Note that the magnitude must be given m weight units, 
not in mass units. 

XVII-XXII Constraint Control Data 

There are 5 types of constrains which can be specified. 

Each constraint type may have different truncation control, 

- although the method used is identical for all types of con- 
straints. The truncation strategy is similar to the one 

used m ACCESS-1, but the sign of feasible region is reversed. 

th , 

If a g constraint function at a design a is evaluated as 

h^ (cT) , h^ (ct) is compared with a truncation boundary value 

(TBV) which is determined by 
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XVII. 


XVIII. 


TBV = + {Mm [h (a) - C] }x TRF + C, 

g q 

where Min xs applied to all g's in the constraint type. 

q 

Initially, TRF is set to be TRF-mitial and at the end of 
each design state, TRF is updated by 

TRF = TRF x (TRF multiplier) 


Since TRF-multiplier is chosen to be less than 1, TBV is 
decreased stage by stage, which means more and more con- 
straints are truncated as the design proceeds. 

Side Co nstraint Control Data 

Since side constraint codes are specified m the element 
data, only truncation control parameters are specified: 


TRF-mitial 
TRF-max 
C-cutof f 


Initial truncation factor 
Upper limit of TRF 
Cutoff base value 
TRF modification multiplier 


TRF-multiplier 
Displacement Constraint Control Data 
Card 1 (15) 

NDPC : Number of constrained displacement degrees 

of freedom 


Card 2 ( E10.4) 

TRF-mitial 
TRF-max 
C-cutoff 
TRF-multiplier 


Initial truncation factor 
Upper limit of TRF 
Cutoff base value 
TRF modification multiplier 
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XIX. 


Min . Norm Ftr . 


Card 3 ~ (NDPC+2 ) 
Node 

1 

Ixyz 


Code 


: Minimum constraint normalization 

factor. Constraints are usually 
normalized by the absolute values of 
the limiting values. 

(315, 5X , 2E10.4) 

: Node number associated with the i-th 

displacement constraint 

: Direction identifier 

0 = not used 

1 = X direction 

2 = Y direction 

3 = Z direction 

: -1 = Lower bound only 

0 = No constraint 

1 = Upper bound only 

2 = Both 


Lower Bound : Lower bound of the displacement 

component 

Upper Bound : Upper bound of the displacement 

component 

Slope/Relative Displacement Constraint Control Data 
This constraint type is restricted to place bounds on relative 
displacement components of two arbitrary nodes. In other 
words, the difference between Y-displacement components of 
the L th and U th nodes may be bounded. But the difference 

r.1_ 

between the Z-displacement of h and X-displacement com- 
ponent of U node cannot be bounded. 
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Card 1 (15) 

: Number of slope/rel-displacement constraints 

(5E10..4) 

Initial truncation factor 
Upper limit of TRF 
Cutoff base value 
TRF modification multiplier 
Minimum constraint normalization factor 
(315, E10.4) 


NSLC 
Card 2 

TRF-initial 
TRF-Max 
C-cutof f 
TRF-multiplier 
Mm. Norm. Ftr 


Card 3- (NSLC+2) 
(L) 


Node 


Node number of the node associated 


with the i 


th 


Node 


(U) 


slope constraint 
.th 


xyz 


Node number of the U node associated 
with the i slope constraint 
Direction and code 


0 

1 

2 

3 

4 

5 

6 


not used 
X direction 

Y direction { 
Z direction 
X direction 

Y direction ■( 
Z direction 


r 


relative 

displacement 

slope 


Upper Bound 
Note : 


Upper bound of the slope/rel. displ. 


1. If I = 1, for example, the constraint function is 

Xy 3 


l - (u 


Node 


' J ' Node ^ ) /Upper Bound £ 0 


- U 


62 




where D„„ is the projection of the distance between 

1 U 

node ^ and node ^ to the Y-Z plane. 

3. If lower bound is to be specified, node^) and 

node ^ should be exchanged to transform it to an 

upper bound constraint. 

XX. Stress/Strain Constraint Data 

Card 1 (1015) 

MTYP' 

Code : Stress/Strain constraint code 

Except for element type 3 

-1 = read stress constrain code element by element 

0 = no stress constraint 

1 = all elements m this element type are constramec 

by lower bounds on compression stress 

2 = all elements m this element type are constramec 

by upper bounds on tensile stress or Von Mises 
combined stress (Element Type 1 or Types 2, 4 ,5, 6 

3 = effectively this implies that both codes 1 and 

2 are applied simultaneously 
For element type 3 

-1 = read strain constraint code element by element 

0 = no strain constraint imposed 

1 = maximum strain envelope criteria imposed on all 

elements 

2 = stress interaction criteria imposed on all eleme: 
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3 


Tsai-Az 2 i criteria imposed on all elements 


Card 2 (7E10.4) 


Initial truncation factor 
Upper limit of TRF 
Cutoff base value 
TRF modification multiplier 
Minimum stress constraint 
normalization factor 
Minimum strain constraint 
normalization factor 
Truss Euler buckling control 
factor 

If TEBCF < 0, TEBCF stands for the specified mean 
radius r of the truss element assuming tubular 
cross section. Stress constraint is 


TRF-initial 

TRF -max 

C-cutoff 

TRF-multiplier 

Mm. Stress Norm Ftr. 

Mm. Strain Norm Ftr. 

TEBCF 


0 2 Max{a= llowable -ir 2 Er/2H 2 } 

If TEBCF > 0, it stands for the thickness to mean 
radius ratio of the truss element (r) assuming 
cylindrical cross section stress constraints 


a > Max {a 


c 

allowable 


-tt 2 EA/[4£ 2 * (|)]} 


If TEBCF = 0 , no Euler buckling constraints are 
considered. 

Card 3 - (1415) 

Stress/strain constraint specification for element type 
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MTYP MTVP 

code. Code <0. If all Code are positive, 

no cards are required. 

MTYP 

For each element type with Code = -1, stress/ 

strain code must be given to all elements sequentially 
starting from element number 1. 

Element stress/strain constraint code. 

Stress code 

-1 : only compression side is bounded 

0 : no constraint 

+1 : only tensile (truss only) or Von Mises 

combined stress is bounded 
+2 : both compressive and tensile stress 

are bounded. 


Strain code 

MTYP 

same as Code specification 

XXI. Natural Frequency Constraint Data 
Card 1 (215) 

NFREQ : number of lowest frequencies to be bounded 

NS PACE : frequency constraint approximation scheme 

0 = first order Taylor series expansion with 

respect to linked reciprocal variables 
(linear in the optimization design space) . 

1 = first order Taylor series expansion with 

respect to linked direct variables 
(nonlinear constrain m the optimization 
design space) 
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2 = second order Taylor series expansion with 
respect to linked direct variables 

Card 2 (7E15.6) 

TRF-imtial : Initial truncation factor 

TRF-max : Upper limit of TRF 

C-cutoff : Cutoff base value 

TRF -multiplier : TRF modification multiplier 

Mm.Norm.Ftr . : Minimum constraint normalization 

factor 

Eig. Conv. : Eigenvalue analysis convergence 

criteria (see note below) 

Acc. Gravity : Acceleration of gravity 

If 0-0, American standard unit is 

2 

assumed and replaced by 38(3.0 m/sec . 
Note: Subspace iteration algorithm is used to obtain 

eigenvalues and eigenvectors. Iteration is judged 
to be converged if the relative differences of all 
eigenvalues are less than Eig. Conv. 

Card 3 (15, 2E10.4) 

Code^ ^ : constraint code 

-1 = lower bound only 

0 — not bounded 

1 = upper bound only 

2 = lower and upper bounds 

~tl*L 

Lower Bound : lower bound on the l frequency 

Upper Bound : upper bound on the l frequency 
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XXII. Not used 


Card 1 : supply one blank card. 

XXIII. NEWSUMT Optimizer Control Card 
Card 1 : 

JPRINT : Optimizer printout control 

standard output = 0 

MAXSTG : Maximum allowable number of stages 

MAXRSF : Maximum number of response surfaces 

per stage; i.e. response factor is 
reduced MAXRSF times before the approxi- 
mate problem is updated. 

MAXODM ; Maximum allowable number of one dimensional 

minimization per response surface 

JSIGNG : sign of feasible region 

1 : feasible region is q^{a)iO 

— 

-1 : feasible region is q (a)kO 

Card 2 and 3 : 

EPSSTG : Stage convergence criterion. 

Overall iteration is judged to be converged 

if both of the following conditions are 

til 

satisfied at the end of the P stage. 

|w p - W p _ 1 |/W p < EPSSTG 

|W p _l " W p- 2 l /W P-l 1 EPSSTG 

EPSODM : Unconstrained minimization convergence 

criterion. Convergence is obtained if 
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RACUT 

STEPMX 


ITP 

Power Fr 
Coefficient 
STEPMX-mul 
STEPMX-L.L. 


the relative values of total function at 
the ends of 3 successive one dimensional 
minimzations are not different by EPSODM. 
Response factor decrease ratio 
Maximum step size at each stage. 

All design variable components are 
constrained by 

STEPMX “ ^ STEPMX. 

Initial transition point for the extended 
penalty function 
specify = 0.5 
specify = 1.0 

Maximum step size modification multiplier 
Lower limit on the STEPMX 
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APPENDIX C 


Examples 

Three examples are given to illustrate input data pre- 
paration for various features of the ACCESS-2 computer program. 

(1) 25 bar space truss 

static constraints only 
mechanical and thermal loads 

(2) 18 element wing box 

static constraints and frequency constraint 
mechanical loads only 
aluminum alloy 

(3) 18 element wing box 

static constraints 

mechanical and thermal loads 

composite material skin and aluminum webs 
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Example 1 


ORIGINAL 

OF POOR QUALITY 



a = 0.276 GPa ( 40,000 psi) 

a a = 23.0 x10“ 6 /K (12.8x10- 6 /°F) 


Loading 








Grid 

Point 

Temperature 
K °F 

P 

X 

N 

lbf 

P 

Y 

N 

lbf 

p 

2 

N lbf 

1 

350 

170 

4448 

1000 

44 480 

10 000 

-22 240 -5000 

2 

350 

170 

4448 


44 480 

10 000 

-22 240 -5000 

3 

311 

100 

2224 

500 




4 

311 

100 






5 

311 

100 






6 

311 

100 

2224 

500 




7 

275 

35 






8 

275 

35 






9 

275 

35 






10 

275 

35 







Fig. C— 1 25 Bar Truss Ref: NASA TND 7965 August 1975 












Example 1 Data Card Image (1) 





Example 1 


Data Card Image (2) 
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1 0 
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0.10.0E-14 


DISPLACEMENT C| 
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. ) — 

I 

0. I 1590E5' 


ONST R/l I NTS! 

__ , 
r_ I 


' 1 N l T v! 


I LOW G 


INE&JL. 


JFL.LQH. 

NMASS 


NSLCL 

T j 


0 o 1 


OOOE 3 


] 

i 1 

[ 

L L ' i 

1 

0. 1 

oooeoj o.s 

r 

[oOOEO 

i ~ 

0 . ljooofco 

1 



sy fi -o 4>Sr 


So 


Sf 


40 Si 



□i>o 


E xample 2 

Truss 
CST 
SSP 

Numbers Within Symbols indicate Element No. 
Numbers at Element Junctions Indicate Node No. 



2 

Lumped Masses : 222.4 Kg m/sec (50 lb) 

at nodes 4 & 6 

Material properties : E = 68.9x10^ Pa, N/m 2 (10x10^ psi) 

p = 2800 Kg/m 3 (0.1 lbs/m. 3 ) 
v = 0.3 

CT U = ( cr L I = 6 . 89xl0 7 Pa , N/m 2 (10000 psi) 


Fig. C-2 Wing Box with Isotropic Panel 











Constraints : 


OF POOR QUALITY 


Side 


Stress 

Displacements 

Frequency 

Loads 


A L = 0.645 cm 2 (0.1 m 2 ) 
t = 0.0508 cm (0.02 m) 
t L = 0.0508 cm -(0.02 m) 
a U = [ a L | = 6.89x10 7 P , N/cm 2 

d 

(10,000 psi) 

U 2 k = l u zkl = 5.08 cm (2.00 m) 
nodes 3 ,4, 5, 6, 7. 
= 25 cycles/sec 
Load Condition 1 

P 1 = 22240 N (5000 lbs) 

P 2 = 0 (N) (0 lbs) 

Load Condition 2 

P-L = O (N) (0 lbs) . 

P 2 = 44480 N (1000 lbs) 
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Example 2 Data Card Image (1) 





ORIGINAL PAGE IS 
OP POOR QUAUTYi 










































Example 3 


Q Truss 

A CST 

« — , 

1 — 1 SSP 

Numbers within Symbols Indicate Element No. 
Numbers at Element Junctxons Indicate Node No. 





Q p C 

Isotropic Material : E = 68.9x10 P. , N/m (10x10 psi) 

’ ci 

*3 3 

Properties p = 2800 K /m (0.1 lbs/m ) 

g 

v = 0.3 

cr U = j a L | = 11. 031 6x10 7 P a , N/m 2 (16000 psi) 
= 23.0xl0 _6 /°K (12.8 x10 -6 /°P) 

Fig. C-3 Wing Box with Composite Panel 
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ORIGINAL PAGE IS 
OF POOR QUALITY 


Composite Material 

: E l = 

1.448X10 11 P , N/m 2 (21xl0 6 psi 

Properties 

Erp ~ 

1.172xl0 10 P , N/m 2 (1.7xl0 6 psi 



4.482xl0 9 P , N/m 2 (0.65xl0 6 ps 


II 

0.21 


p 

1549 Kg/m 3 (0.056 lbs/m 3 ) 


“l = 

-3.779xl0“ 7 /°K (-. 21 x10~ 6 /°F) 


°T = 

2 . 8797xl0 5 /°K (16x10 _6 /°F) 


t’ 



% “ 

0.008571 m/m (0.008571 m/in) 


tf* 

o 

II 

-0.008571 m/m (0.008571 m/m) 


t , 

^ ~ 

0.004706 m/m (0.004706 m/m) 


Hf 1 

II 

-0.017647 m/m (-0.017646 m/m) 



0.018462 (0.018462) 

Constraints : Side 

: A L = 

0.645 cm 2 (0.1 m 2 ) 


,L 

t ) = 

0.0508 cm (0.02 m) 


L 

T — 

0.0508 cm (0.02 m) 


Stress and 

strain : A^ given m the Material Properties 

Displacement: = |u^ k ! = 5.08 cm (2.00 in) 

at nodes 3 , 4 , 5 , 6 and 7 


Loads : 

Mechanical Loads : Load Condition 1 

= 22240 N (5000 lbs) 
P 2 = O(N) (0 lbs) 

Load Condition 2 


P-L = O (N) (0 lbs) 

P 2 = 44480 N (10000 lbs) 
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Load Condition 3 


P 1 = -22240 N (5000 lbs) 
P 2 = 0{N) (0 lbs) 

Load Condition 4 

P-j_ = O(N) (0 lbs) 

P 2 = -44480 N (10000 lbs) 

Thermal Loads: All load conditions 311°K (100°P) 
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Example 3 Data Card Image (1) 
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Example 3 Data Card Image (3) 


VO 

H 


15] 

1 

1 


l 6f 


3 


20 

1 

3 

T?| 

3! 

1 

l 


1' 

0. 1000E3 
— a ol 


1 

1 

I 

.. X 

1 

1 

l 

-il 

1 

1 

1 

-? 

I 


?! 

4 

_-4i 

18i 

4 

li 


I* 

l 1 

4 

1 

1 

1 

-1 

1 

1 

1 

1 

1 

1 


2 
4] 

6! 

ao.oo-E-o 


0. iJoOOE O’ 
3| 


1 


0 . 5j0 0 0 E 2 
0.5000E2 
._O-.60l0OEP 

0.6000E0 
2 


5 

6. 

1 


I 


3 1 
0. X 
0] 


o 1 

jOOOE 0, 
0| 


0] 

0 • 0 

1 .0 


J0- 

l 



SoiOEO] 

DOOEOi 


0 .doi OEOj 
O.X200E1 1 


19 

51 

1 

1 

1 

1 

1 

- I 
1 
X 
1 

. X 
1] 
X 
X 
X 
l 


— &!_ 

20 

6! 


__X_. OP 0 QF-0| _ 1L. 2SQ£L0H0i. . 

i i i 

|00 OE It 1.2|OOOEO| 


0.1 

-0 .2 


00 0 El 



it - 

li 

1 

-l - 

i. 


1 


___z 

1 

! 

Hi.. <3 

_ J_CL 

. lT 

i 2 

13 

r 

1 

1 * 

li 

7 

8 

li 1 

1 1 ] 

| 

, 1 . l! 

_ L 

1 

_ l 

! 

X 

ii ._ 1 

1 

1 _ . _ L 

n 

1 

__ij 

i 

i 

ll 

L 1 

l i 

.. li 
1 

1 

l 

( 


i 

_i 

l 

l 

, . 1 

>i 

li 

1; 

_ ii 1 

1 


, - 1 

1 


r 


L 

-t- ■ 


1 




- 0.2 

000E1 

- 0.2 

000 El 

- 0.2 

00 OE 1 

! 

0 

aj 

1.0 

OOOEO 


pan of if 

-0^2|0.OP -JXm ^ 00()El | 

0.2100 0 El, 
00 OEX, 


| 0.200 1 

r -) - 


1 • 0 00 0E0 




1® 2,00 0E0 0.1 000E5; 


1 4 


i - 


_l L 


I FLUSH,. 


NSLC 


oooeo! o.io 


NQ FREQUENCY. CONSTRAINTS 


OE-14 Oo 1 00 0E3' O.lOOOEOi 0.5000EO, O.l.OOOEO 







30 


3 <~ 


YO 


y$ 


£0 


7<r 


3o 



APPENDIX D. Output for 25 Bar Space Truss Example 
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RASIS CUTOFF FACTOR 1 .0000 
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STEP SIZE MODIFICATION FACTOR 0.10000*01 
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0.291 020 +00 


26 SIDE CONSTRAINTS 
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O.R0670+00 0. 

1 0930 +0 1 
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t- 1 DISPLACEMENT : -l^LQWEE BOUND. +1=UPP C R SOUND 

O SLOPE i ALWAYS UPPER BOUNDS 

|_i STRESS/5TP A I N * I O-VON WISES 

i=longttudinal strain 

2=TRANSVERSE STRAIN 
3=SHCAR STRAIN 

FREQUENCY ASSOC IATFO MODF NUMBER O® 




ANALYSIS TIME DATA 

ASSEMBLE MASS/STIPFNFSS MATRIX 

ASSFMBLE LDAP VECTORS 

DFCOMPOSE STIFFNESS MATRIX 
SOLUTION OF DISPLACEMENTS 
FREQUENCY ANALYSIS 
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n - ACAgAAn-A I 
0 . 1 0520BD-0I 
0 .9 I6667D-02 
0 .0 


CONST PAINT EVALUATION 

POSTURE TABLE SET 

SELECTIVE GRADIENT EVALUATION 


0 . 7281250-0 1 
0.645833D-0I 
0 .S73A37DFQ0 
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CONTROL PAR \ MET rrs 


INITIAL transiticn pcint r.n = ,1.13'iont^ 

T F AN S r T I ON POINT EXPONENT P = 0 . = 3 C C'' + 00 

INITIAL TRANSITION POINT CIEFFICIfnT C = 0.19000+90 

GOLDEN SECTION CONVERGENCE gPSGSN = 0.13 n ''r-32 

UNCONS TRAINED M ININ I TATI 3N CONVERGE NC r EPSODM = n , I poor- 3? 

CONVERGENCE among RFSDJN5E SURFACES SPCRSF = 0.5JCCC-92 

RESPONSE FACTOR PFDUCTION RATIO NACIJT = 0.300CF+00 

MINIMUM ALl CWABLE RESPONSE FACT3R RAMIN r D.lDCOP-IA 

M _A. X I NU M_A LL C.W A B L g._ S T EP , .5 ,Jgg ST QP_MX — E — 0, 1 TOO OFT j 

MAXIMUM ALLOWARLF GOLD 1 " N SECTIONS MAXGSN * 20 

MAXIMUM NUMPFR CD 0.0. M. P r R SURFACE MAX0DM = E 

MAXIMUM ALLOWABLE RF S°ONSE SURFACES MAXP6 + = 2 

PF1NTOUT CONTROL .IF.INT = O 

FII-ITE DIFFFPCNCF GRADIENT CON T ROL IFt = C 


SYSTEM PAPAMl T= RS 

NUMBER OF DFSIGN VARIABLES NDV = ?F 

NUMPFR CF EFFECTIVE CONSTRAINTS NTCt = 22 

INITIAL DESIGN ANALYSIS SUMMARY 

IN IT I AL DESIGN VARIABLE VECTOP 



0.1 0000*01 
0 . 10000+01 
0 . 10000+01 

0 . 1 0000 + 01 
0 . 1 0000+01 
- 0 . 1 0000+01 

0 . 10000+01 
C. 1 0000 + 01 
0 .10000 + C 1 

0 . 10900+01 
C, 1 0000+01 
0. 1 00 0 0 + 0.1 

0 . 1000 D +01 
0 . 1 0000+0 1 
0 . 1 .C 00 D +0 1 

0 . lOACQ+O 1 
0 . 10000+01 

0 . 10 D 0 D +01 

C. 10 C 90 +D 1 

0 . 1 0000+01 
0 . 1 009 D+ 0 1 

C. 1 D 03 D+C 1 

O.lCPOC+Ol 

3 » i ^>o*n+n i 
0.1 0 C 0 rx + 01 


S IDF CONSTRAINTS 










-I 

1 

-2 

2 

—3 

3 

-4 

a 

— 

s 


0 *99000 *00 

0 . 9 R 00 r »+02 

0 . 99000*00 

0 . 9 O 90 D +02 

0 . 99000+00 

0 . 99000+02 

0 .$$000 +00 

0 . 90000 + 0 ? 

o .'9 "0 jr+ rn 

0 »ooooo*r? 


-6 

6 

-7 

7 

-6 

6 

—a 

9 

-.10 

t 0 


o.gscoc+oo 

O.RPOOO +02 

0 .99000 *C 0 

0 » 9 Q 30 D + 02 

0 • 9 V 00 P *00 

0 . c Q C 0 P + 0 > 

0 .' 3 OC? 0+*0 

0 .CC 0 OP+C 2 

C. K'C ^1 *-o ) 

D *oo 000+0 2 


-i l 

1 t 

-l? 

12 

~t * 

! 3 

- I'i 


-l * 

1 * 


0 . 99000+00 

0 .RO 0 CP +02 

C • <90000 + 0 0 

P . 90 90 0 + 02 

9 • ^ 9 non *00 

O. 9 ^ 900 + J 2 

Of JOC ^D+DO 

D .OooOD+O? 

0 .O 9 CCD+ CO 

0 .^^ 000+^2 


-16 

16 

-17 

l 7 

-l* 

l S 

- l« 

1 9 

-2 y 

?0 


O. 9 C 0 OO *-00 

0 . CROC D +02 

0 . 9 = 000+00 

0 . 99000 + 0 ? 

0 .uaoon *°0 

c. gcooo + 'V 

o . 30 +oo 

r . 99000 + 6 ? 

C. O 3 CCf' + ^0 

0 .<3 ^ CO + C ? 

i — 1 

-21 

21 

-22 

22 

“23 

23 

-?* 

24 


?S 

o 

C.QQOOO+CO 

CONSTRAINTS 

o. Tooor+o 2 

0 .90000 + 00 

0 . 9900 D +02 

0 . O 900 P +00 

0 . <= 9990+92 

0 .ocrjnn 4-00 

O.o-OOD+O? 

0 .o< 0 0 L+ 00 

0 . 00000+02 


0 . e 750 O+C 0 

0 . 30000 + CO 

C. 9429 C+ 0 C 

0 . 9109^*00 

o . ppeco+cc 

0 . 90970**0 

0 .rtQ 190 * 0 ^ 

C • 2 1 70+0 0 

C. - ^0 *IT* CC 

0 • A 4 ./J 30 + 00 


0 . 961 1 0+00 

0 . 93960 + 00 

0 . 96050+00 

0 . > 3330+00 

0 . 90670+00 

0 « $ 0440+0 > 

C .'JbA ju+no 

0 . 9 ^ 5 T P +00 

o.^^arL + oo 

0 . 40 + 00 


<3.e704D+00 O.S42AP+0O 

OB JECTIVE FUNCTION = 0.656 13270+03 

DIRECTION FINDING 

TRANSITION POINT = 0. 1 000000 + 00 COEFFICIENT PF TRANSITION POINT CALCULATION = 9.27447AC-91 


DIPECTICN finoing 

TRANSITION POINT = 0.1662640-01 COEFFICIENT OF TRANSITION POINT CALCULATION 

PIRFCTION FINCINC 

TRANSITION POINT = 0. 1972960-01 COEFFICIENT or TRANSITION POINT CALCULATION 

0 I RFC T ION FINDING 

TRANSITION POINT = 0.201905D-01 COEFFICIENT OF TRANSITION Ft" I NT CA1CUI.AT10M 

DIRECTION FINDING 

TRANSITICN POINT - 0 . 2030 24 C— 91 COEFFICIENT nr Tf ANSI r IOM POINT ' AL CUL A T I ON 

===== D! RFCTTCTR FINDING ~ — — — — 

TRANSITION POINT = O.Y 100230-02 C 1EFF I C I FNT JF TRANSITION FP I NT CALCULATION 


O.c.t 1 ?47r-02 
0 .^4 l 6270-02 
O -6Fa ) 7pn-C7 
' .'■•7?sor-n? 
0 .iogaoap-02 


! 


DIRECTION FINDINC 
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FINAL RESULTS CF OPT I M I Z A T I CM 



0.943CD+00 0»6A3<iD*00 0.SS89P+00 0.^7440+00 0.67093+00 0.603^0 + 00 0.74T3D+C0 0 . 741 TD + C 0 O.ASiOO+OE 0 « 3 7-) <ir> 

C.46C4O+C0 0.367SD+C0 


TOTAL FUNCTION = 0. 284004 10+03 

no i cr t n/c pi imtt t om — fi . i t 


FINAL statistics 


NUMBER OF ONF DIMENSIONAL SEARCH . 

number of analyses 


approximate OnjECTIVF FUNCTION ...... 

GP AD IF NT OF OBJECTIVE FUNCTION 

CONSTRAINT FUNCTIONS 

GRAOrr.NT OF L f N F AQ CONSTRAINT FUNCTIONS ... . 
GRADIENT OF N0NLUE AS CONSTRAINT FUUCTI JNS 
APPROXIMATE CONSTRAINT FUNCTIONS 


C . 0 * M . 
CIRFCTS0N 
EQ. SOLVIT 


CUMULATIVE CPU(<=EC1 


0,2*04 

O.SilR 
0.3' 17 

rt . A r 


0,0 

0 . 00?7 

o . o i o s 
0-0001- 
0.7 

n .02 7P 







STAGE NO 


2 


APPROXIMATE PROBLEM generator 


CURRENT MEMBER SIZE 


MEMBER TYPE NUMBER I 

0.46380-01 0.50580 + 00 0.48930+00 0. 10220 + 00 0.40120 + 00 0.4+670+00 3.M’)T'+0'' 0.4*700+00 0.37*40+00 0.57*40-01 

0.I278D + 00 0.29010 + 00 0.14920+00 P. 20120 + 00 0. 131 70 + 0 3 0.23380+00 0.508*0-91 O.3704r+OO C.37A81+9C 0.21010+00 

0.2041D + 00 0.59510 + 00 ,0 . 6.6.110 + 00 0.51590 + 00 Q_. 7 £3 2 0+0X1 


CURRENT WEIGHT DATA 


MFMBER TYPE NUMBER 1 W F I GMT = C. 1 131500+03 


VARIABLE STRUCTURAL WEIGHT 0.M3150D+03 

FIXED STRtlCTURAI WF I.XIHI 0^0 

TOTAL STRUCTURAL WEIGHT 0.1131500+03 

NON-STRUCTURAL WEIGhTHT 0.0 


THTAI WEIGHT O. 1 I 31500+01 

CONVERGENCE CHECK STAGF NC.= 2 0.4B16D+01 0.151Rp+2e MUST tt c LLSS THAN 0.100000D-02 

OBJECTIVE FUNCTION OF THREE CONSECUTIVE STAGES ARF O.lOOOOCD+31 0.6S31T1P + 03 0.1131500+01 


POSTURE TABLE 



R«T A IM r D 

TOTAL MPMREP 

MEMBER 

NOr r QIF^rTjQfs 

L . C . 

MOHF 

C CN 3 T P A IMT VALUl S 


<* 

SIDE CONSTRAINTS MOST CRITICAL 

= 0. 784401 0+00 




M 

1 

i 

i 

1 



-10 

0* 754 \n\ I'M 00 


O 

2 

4 

i 

4 



- 1 0 

0, '*66™ 7U + 0 0 


cn 

3 

Q 

i 

J 



- I 0 

0. 973432^+^0 



4 

1 0 

1 

i n 



-in 

— n.pqio4 



5 

1 1 

i 

1 1 



-in 

0. 021 ’330+00 



6 

1 2 

i 

1 2 



-10 

C.Qf6 { iA(| P+00 



7 

1 3 

i 

1 3 



-10 

0.9329*20+00 



B 

1 4 

1 

1 4 



- 1 « 

MKKMUELnXiH 

(v- 


Q 

1 s 

i 

1 5 



-10 

0.924 04 3D + 0 0 



10 

16 

i 

1 6 



-10 

0. 957232C+0C 

% s 


l i 

1 7 

1 

1 7 



-10 

0 . *03 J7ID+00 


i p 

t-P. 

« 

\ P 



....... - 1 o 


O t3 


13 

19 

i 

1 R 



-10 

O.973JO1L+0C 



14 

20 

i 

20 



-10 

0.9523-60+00 

td > 


15 

2 1 

i 

2 t 



-10 

0. Ccoqn^MOO 

tr* 
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CONSTPAI NTS 

MOST 

CRITICAL = 

0 . 4* 77* PQ + C C 

1 





16 

27 

1 

1 


1 

0 

0 .+F4461U+00 



17 

28 

I 

2 


1 

n 

r. c :>6»*PO+00 

ft 


tft 

TO 

1 

-* 


1 - 

r> 



19 

J 3 

1 

4 


1 

0 

C . + *2 7<+ »D+«0 



20 

35 

1 

5 


1 

0 

0* C 2<?0^4p + n 3 

l 3 t-j 


21 

36 

1 

6 


1 

n 

0 .5308*30+00 

ts CO 


22 

39 _ .. 

1 

7 


L. 

n 

n . ^4. SQ7 AI 0 .+. no 



23 

40 

1 

S 


i 

0 

o • c a 20 7 jr+ oo 



24 

43 

1 

9 


1 

0 

0 « f 00 f Oc 



25 

45 

1 

1 0 


i 

0 

0 . °S7 ?1hD + 0C 



26 

47 

1 

1 1 


i 

0 

0. 84°094t +00 



















4. 

1 

J l •* 

1 

" ♦ ■ ♦ / 

23 

55 

J 1 5 

1 0 

0*^41 t <V3 

29 

*6 

1 l f3 

1 0 

0.‘ 1C4&U)* -> 


ri_g 

1 17. 

1 0 

n - 7! 

31 

eo 

I 1 3 

1 r 

c ♦ T i d<ac^r<-c 

32 

<=.2 

J i 9 

t p 

c . f OC7 SOP4 c 

13 

65 

J 20 

1 0 

C . 7 l‘-A,Vin4 ^ 

34 

67 

! 2 1 


mmm wan mi 

35 

69 

1 22 

l c 

C • * 7 7/ 0°f >♦ 0 

36 

70 

1 23 

1 0 

0 . 4 A 777 JfHQ 

37 

73 

1 24 

t 0 

0. < 7221-204 0 
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NODAL PISPLACFMPN IS 


NnnF x X. 


Z 


ilQOB. 


X. 


JL 


LOAD CONDITION 1 


1 -0 . 176290-01 0 . 35760C + 0I 0 . 1 3 0 7 2P * 0 t 

3 0.138370*00 -0.233960+00 -0 . 78 28 3D* 00 

S 0 .221020+00 -0 .3031 70 + 00 0.122850 + 01 

2 Q-*£ QjlQ OjlO 

9 0.0 0.0 0.0 


2 -C.17A72r-"l I.Jinil tni -''.-(.2031+00 

4 -A.-(?i?or-i)i - o .2 ja i 7 r +r o -*■ . ’<-f '■'o+oo 

t ' -0.170960+00 — 3 *3 fi 3041 +0 C 0.1197(11*01 

_a o . n !LlC i^o 

10 c . o o.o o.o 


25 SIDE CONSTRAINTS V TO 25 MOST CRITICAL CONSTRAINT^ 0 . 1 3 7*= 9530-01 

0. 18030-01 0 .95430 + 00 0.95660+00 0.eS49D+00 0.9467D+00 0,9=890+00 0.95050 + 00 0.95430+00 0.93T1C+00 0.15160-91 

0 . 15060-0 1 0.36 930+00 0...8121D + 00 0.69BQD+00 8.84340+00 0. 90 33.D_+_Q0 0.7.1930 + QC 0.66S3D.+.00 0.7314Q + CC-. 0. 1.1760-01 

0.218ID + 00 0.971 0D + 00 0.9766D+00 O.S66CD + 00 0.97750 + 00 



WTYP 

M 

LC 

s-cnMniN^n 

sx 


SY 

Sky 

^Y+-THf P . <;y-THCrv RXY-Tl r C.M_ 














1 

l 

1 

-0. 430025D+0A 

-0 .43000+04 

0.0 

0.0 


- 



I 

1 

1 

-0. 4300250+0* 

-0.43 00 0+0 4 

0.0 

0.0 





1 

2 

1 

_ -0 . 3O90P0D+06 

-0. 39 910 + 0*=; 

n.n 

n. o 





l 

2 

1 

-0. 3O9050D+05 

-0.3991D+05 

0*0 

0.0 





1 

3 

1 

-0. 3987S1D+05 

-0.398*0+05 

0 . 0 

0.0 





1 

3 

1 

-0. 39 87 S ID + 0 15 

-0.39880+05 

0*0 

0.0 





1 

A 

1 

0.3985420+05 

0 .. 19359 + 05 

n _ o 

n - f> 





1 

4 

1 

0 . 3985520+05 

0 *39 6 50 f 05 

0.0 

0.0 





1 

5 

1 

0. 39 86 7 1 D + 05 

0 * J907D + O5 

0*0 

0.0 





I 

5 

1 

0. 3O8671D+05 

0 .39870 +05 

0.0 

0.0 





1 

6 

1 

-0. 39896 90 *05 

-0.39900+05 

0.0 

0.0 





1 

6 

1 

-0. 3089680+05 

- 0 • 39900 + 05 

0.0 

0.0 




1— 1 

1 

7 

1 

0. 39 89 7 OD + 05 

0 .39900+05 

0.0 

0.0 




o 

1 

7 

1 

0.3969700+05 

0.3990D+05 

0.0 

0.0 




CO 

1 

3 

1 

-0. 3986960+05 

-0.39870+85 

-.0 .0 

0.0. 





1 

8 

1 

-0. 3986 060+05 

-0.39970+05 

0.0 

0.0 





1 

9 

1 

0.3963560+05 

0.39840+05 

0.0 

0.0 





1 

9 

1 

0.39835*0*05 

0 .39840 +05 

0.0 

0 .0 





1 

10 

1 

0. 4411 520 + 04 

0.44 1 20+04 

0.0 

0 .0- . 





1 

1 0 

1 

0*4411 520+04 

0.44120+04 

0.0 

0.0 





i 

t I 

1 

0.4396MD+04 

0 .43990 +04 

0.0 

0.0 





i 

1 1 

1 

0. 4398910+04 

0 .4 39 40 + 04 

0.0 

0.0 





i_ 

12 

1 

-0. 3951920+05 

-0 .39520+05 

0-0 

0.0 





i 

12 

1 

-0.3951 920+05 

-0.39520+05 

0.0 

0.0 





i 

1 3 

\ 

0. 391998D+05 

0 .3920D+05 

0.0 

0.0 





i 

1 3 

I 

0.39190PD+05 

0.39200+05 

0.0 

0.0 





t 

1 4 

... 1 

-0. 39e563D+05 

-0.39860+05 

0 . 0 

8.0 





i 

14 

1 

-0.3985630+05 

- 0.39860 + 0? 

0.0 

0.0 





l 

1 5 

1 

0 - 3960330 + 05 

0 *3^JR OD *05 

0.0 

0.0 





i 

15 

1 

0. 3980330+05 

0 . 39800+05 

0.0 

0.0 





i 

16 . 

.. 1 

-0. 3987540+05 

-0 .39080+95 

0.0 

0.0 





i 

16 

1 

-0 . 3987940+05 

-0 • 39030+05 

0.0 

0.0 





i 

1 7 

1 

0. 3960670+05 

0.39810+05 

0.0 

0.0 





i 

17 

1 

0* I960 570 + 05 

0.39810+05 

0.0 

0.0 





i 

l 8 

1 

-0. 2999570+05 

-0 .3000D+05 

0.0 

0.0 





I 

18 

1 

-0. 299^70 + 0^ 

-0.30000+05 

0.0 

0.0 





i 

1 9 

1 

-0.33A2C7D+05 

-0 .33420+05 

0.0 

0.0 





i 

19 

1 

-0 . 3342 07D + 05 

-0.33420+05 

0.0 

0.0 





1 

20 

1 

0. 31 144 70 + 05 

0.31 140+0? 

0.0 

0.0 






1 

<■ u 

i 

• J 1 1 " ' ' * *7 


t 4 


1 

21 

1 

0. 27=7430+06 

0 .27 )70 + 0 S 

0.0 

0.0 

1 

21 

1 

0. £197 A 30 + 05 

0 ,27970 +05 

0. " 

0.0 

I 

2? 

1 

n . iQP. 

0 - -np 4-n r*. 

_ja^a 

M ^ 

1 

Z2 

1 

0 • J <5034 7D + 0* 

0 .J ?<3 JO *-05 

0.0 

^ « 0 

1 

23 

1 

-0.3°B74 70 + 0 = 

-0 , J9R7D *-05 

0.0 

0.3 

1 

23 

1 

-0. J = 8"A 70 + 05 

-0 .39F370 +05 

0.0 

0 , A 

1 

-2A 

_L_ . 

n- 3££31i?D+f>5 

n „ ?Q.q »• 0 5 

0 - 0 

■>- r> 

1 

24 

1 

0. J9P31OD+06 

O.T«8JD+OF 

0.0 

0.0 

I 

25 

1 

-0 . 3909300+05 

-0 .30040 + 05 

Ot 0 

0*0 

1 

25 

1 

-0. 3 9 80 300 + 0 s 

-0 . J9890+OS 

0 • 0 

0.0 


50 STPCSS/STRAIN CONSTRAINT? 26 TO 7S (OST CrlTTCAL CTOiTPAIMT; 3 . J376 67 70 - ” 2 


0. 8925D+00 
0. P580D-02 
0. 1 1 1.0D + 0 1 

0.11 080+01 
0.10070+01 
0 .80008 + 00 

C #23750-02 
0 . 1 $070+01 
n . 1 P02D-0 ! 

0* I9QR0+0 1 
0,25750-02 
n- i_g.R.qn + r>i 

0.31240-02 
0 .32590-02 
n . 1 0400+01 

0. 19970 +01 
0. 1 9970 *01 
n. znotr -01 

0 ♦ 1 0 + C 1 

0* l c ^r.f 0 1 
n * iboba- a? 

0. 

i 1 id -02 

O , 1 MAr 1- n 1. 

9 . 1 9C7Of0 1 
c.mouM 

n^; Q£f.i>* r. t 

0 *3 V4p- 3? 
0 *05 ?7C + 00 
r* . ^ 0 1 r i r>~ rv j 

0.301 oo-oz 
0. If O9D + 0 1 

0 .19470 + 01 
0.30 06 t> + 0C 

0 .19950+01 
0. 1=960+01 

0.476TD-02 

0.41320-02 

O.ZbOlO+OO 
0 . 31 J20-0? 

0. 1 7500 *01 
0.19470+01 

0 . it 4&0+00 
0 « l OQ 6 ( ) + 0 1 

0* 1 * ;J6Cf 01 

o.A^^br-02 

Q« j 7?<'r, + Q J 
0# 26750-02 

6 *221 4n+uo 
0 .19970+01 


■CCNSJ-HAINTS OUT l~)F 2 S CIITrFF POINT:; fl . 4nR?S6n+nn 

1 to 11 20 21 


0 CONSTRAINTS OUT or 0 CUTOFF P01NT= 0.40SPS6C + 00 


0 CONSTRAINTS OUT OF 0 CUTCFF PO I NT - 0.4082*60+00 

2? CONSTRAINTS OUT OF e 0 CUTCFF POINT* 0 . 40 1 4250 + 00 

ZB an 3J as 3£ Jfi jlC _++J _AH 34- 

52 5*= So 59 PO C2 OS 57 C9 70 

7 3 7 A 

?? CnNSTRATNTS OUT OF ~PP PFTATNFD nilF TO V AP X AFt F -LI MK 1 N O 

28 30 33 Jo 36 JO AO 4J ho SI 

S2 5S 86 5? 60 62 65 6" *0 70 

73 7 A 


0 CONSTRAINTS OUT or O CUTCFF POINT* 0. 4014250+00 

0 CONSTRAINTS OUT OF 0 CUTCFF POINT* 0,4014250+00 


****************************************** A**** *** INTEGER ARRAY" SI7r 20000 .G r , 
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